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ABSTRACT 

The  production  of  nitric  oxide  by  soft  solar  X-rays  in 
the  E-region  of  the  ionosphere  is  investigated.  An 
empirical  expression  for  the  variation  in  X-ray  flux  as  a 
function  of  F10.7  is  determined.  This  expression  is 
incorporated  into  a  one  dimensional  diffusive  photochemical 
model  to  compute  nitric  oxide  densities.  The  results  of 
these  calculations  are  compared  with  NO  observations  from 
the  Solar  Mesosphere  Explorer  satellite.  Variations  of 
X-ray  flux  by  a  factor  of  30  over  the  solar  cycle  can 
explain  the  observed  variation  in  nitric  oxide  densities. 
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INTRODUCTION 


I . 

This  thesis  investigates  the  production  of  nitric  oxide 
(NO)  by  solar  soft  X-rays  in  the  E-region  of  the  ionosphere 
over  the  solar  cycle.  The  variation  between  extreme 
ultraviolet  (EUV)  flux  and  X-ray  flux  as  a  function  of  solar 
cycle  is  investigated  as  the  principal  factor  influencing  NO 
production  in  the  E-region.  This  research  was  prompted  by 
the  hypothesis  of  Barth  et  al.  (1988)  that  nitric  oxide 
variations  are  produced  by  variations  in  the  solar  soft 
X-ray  flux  between  20  and  100  A.  The  main  topic  will  be  an 
analysis  of  the  production  of  nitric  oxide  at  110  km  as 
measured  by  the  Solar  Mesosphere  Explorer  (SME)  satellite. 
These  measurements  are  compared  to  the  results  of  a 
one-dimensional  diffusive  photochemical  model  developed  by 
Cleary  (1985)  at  the  University  of  Colorado  and  modified  by 
Siskind  (1988)  to  allow  for  large  variations  in  the  solar 
X-ray  flux  between  18.62  and  50  A. 

A.  THESIS  GOALS 

This  thesis  is  a  theoretical  and  observational  study  of 
the  production  of  nitric  oxide  in  the  E-region  of  the 
ionosphere  based  on  variations  between  EUV  and  X-ray  flux 
over  solar  cycle.  An  empirical  expression  relating  X-ray 
and  EUV  flux  is  sought  in  order  to  better  model  the 
production  and  loss  of  NO  in  the  lower  thermosphere. 
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B.  THESIS  OUTLINE 

The  thesis  is  divided  into  seven  chapters.  Chapter  II 
presents  background  information  about  the  composition  and 
structure  of  the  ionosphere,  and  describes  some  techniques 
for  measuring  the  density  profile  of  various  constituents. 
Atmospheric  chemistry  is  reviewed  in  Chapter  III,  with  the 
emphasis  on  photochemistry.  The  continuity  equation  and 
physical  transport  mechanisms  are  discussed,  as  well  as  the 
method  used  by  the  model  in  solving  the  continuity  equation. 
The  chemistry  of  nitric  oxide  is  focused  on  in  section  C, 
and  the  concept  of  X-ray  scaling  is  introduced.  Chapter  IV 
discusses  photochemical  modeling  and  describes  the  basic 
structure  and  operation  of  the  photochemical  model  used  in 
this  thesis.  Chapter  V  proposes  an  empirical  relationship 
between  the  X-ray  scaling  factor  (FRACXR)  and  10.7  cm  flux 
over  solar  cycle,  with  results  of  comparisons  between  model 
data  and  SKE  satellite  measurements  presented  in  Chapter  VI. 
Chapter  VII  summarizes  the  thesis  and  presents  suggestions 
for  further  research. 
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I I .  BACKGROUND 


Nitric  oxide  makes  up  only  0.5  parts  per  billion  of  the 
Earth's  total  atmosphere  ( NOAA ,  1976),  yet  its  role  as  a 
catalyst  in  the  destruction  of  stratospheric  ozone  makes  it 
one  of  the  most  important  trace  constituents.  The 
ionization  of  NO  by  geocoronal  Lyman-a  is  largely 
responsible  for  maintaining  the  nighttime  E-region 
ionization  at  about  100  km.  Additionally,  because  nitric 
oxide  is  very  reactive,  it  plays  an  important  role  in  the 
daytime  ionospheric  character.  Because  this  thesis  concen¬ 
trates  on  ionospheric  NO  production,  a  brief  description  of 
the  ionosphere  is  presented. 

A.  THE  IONOSPHERE 

Above  about  60  km,  there  are  sufficient  ions  and 
electrons  to  create  a  region  which  affects  the  propagation 
of  radiowaves.  This  region  is  called  the  ionosphere .  The 
ionosphere  is  neutral  overall.  It  is  composed  mostly  of 
neutral  molecular  and  atomic  gasses  which  have  concen¬ 
trations  that  are  several  orders  of  magnitude  greater  than 
the  ions  and  electrons.  The  variation  of  electron  density 
with  altitude  leads  to  three  natural  subdivisions  within  the 
ionosphere.  These  three  regions  are  called  the  D,  E,  and  F 
regions.  The  F-region  is  further  divided  into  two 
subregions  called  the  Fl-  and  F2-regions.  This  character- 
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Altitude  (km) 


istic  "layering"  of  the  ionosphere  is  due  to  the  variation 
of  atmospheric  composition  with  altitude  and  the  wide  range 
of  ionizing  wavelengths  associated  with  each  neutral 
species.  Figure  2.1  shows  the  different  regions  of  the 
ionosphere . 
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Figure  2.1.  Regions  of  the  Ionosphere 


1.  The  D-Region 

The  D-region  is  the  lowest  lying  ionospheric  region 
and  is  therefore  produced  by  ionizing  radiation  at  those 
frequencies  which  penetrate  deepest  into  the  atmosphere.  It 
extends  from  about  50  to  90  km  in  altitude.  The  major 
source  of  ionization  in  the  D-region  is  photoionization  of 
NO  by  the  Lyman-a  1215.7  A  line. 

2.  The  E-Region 

The  E-region  lies  between  about  90  and  130  km,  and 
has  a  daytime  electron  density  about  two  orders  of  magnitude 
greater  than  that  of  the  D-region.  Ionization  in  this 
region  is  primarily  due  to  photoionization  of  molecular 
oxygen  and  nitrogen  by  EUV  and  soft  X-rays. 

3.  The  F-Region 

In  the  F-region,  above  130  km,  O,  C2,  and  N2  all 

photoionize,  with  the  primary  source  above  150  km  being 

atomic  oxygen.  Within  the  region,  the  electron  density  is 

fairly  uniform  until  about  160  km,  where  it  then  rises  to  a 

6  3 

daytime  maximum  of  approximately  10  /cm  at  220-450  km. 

Because  of  this  characteristic  profile,  the  F-region  is 

further  divided  into  the  Fl-  and  F2-regions  lying  at  130-160 
km  and  160-600  km  respectively.  In  the  Fl-region,  ion 

chemistry  dominates  the  characterization  of  the  ionosphere 
whereas  in  the  F2-region  physical  transport  (i.e.,  dif¬ 

fusion)  is  the  most  important  process. 
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B .  THE  THERMOSPHERE 

The  Earth's  neutral  atmosphere  is  divided  into  regions 
primarily  on  the  basis  of  temperature  gradients.  Typically, 
at  middle  latitudes,  the  temperature  falls  with  increasing 
altitude  until  reaching  a  local  minimum  of  around  200  K  at 
about  10  km.  Here  the  atmosphere  undergoes  a  temperature 
inversion  where  the  temperature  rises  with  increasing 
altitude  until  a  local  maximum  of  290  K  is  reached  at  about 
50  km.  The  atmospheric  temperature  then  again  falls  with 
increasing  altitude  until  an  absolute  minimum  of  about  160  K 
is  reached  at  90  km.  This  minimum  is  termed  the  mesopause 
and  marks  the  beginning  of  the  thermosphere .  In  the  thermo¬ 
sphere,  temperature  again  increases  with  altitude, 
approaching  a  maximum  value  in  the  region  of  1000-2000  K  at 
a  height  of  several  hundred  kilometers.  Above  the 
thermosphere,  the  atmosphere  is  essentially  isothermal. 

C.  REMOTE  SENSING  TECHNIQUES 

Because  the  ionosphere  is  a  plasma,  it  will  affect  the 
propagation  of  radio  waves.  It  is  for  this  reason  that 
radio  sounding  methods  were  the  first  experimental 
techniques  used  for  studying  the  ionosphere  (see  e.g., 
Bauer,  1973).  Other  techniques  can  be  used  to  obtain 
experimental  information  about  the  chemistry  and  chemical 
processes  of  the  atmosphere.  Among  the  most  important  of 
these  are  ground-based  observations  of  airglow  emissions  and 
the  use  of  satellite- ^orne  photometers  and  spectrometers. 
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1.  Remote  Sensing  of  Chemical  Processes 

a.  Ground-based  Airglow  Observations 

The  term  airglow  is  used  to  describe  non-auroral 
light  emissions  from  the  upper  atmosphere.  These  emissions 
are  extremely  weak  and  therefore  require  very  sensitive 
equipment  to  be  detected.  A  common  instrument  used  for 
airglow  measurements  is  the  scanning  monochromator  which 
detects  light  in  a  very  narrow  bandwidth.  A  major 
difficulty  encountered  in  airglow  detection  is  the  isolation 
of  the  emission  from  the  strong  background  of  scattered 
sunlight,  (see  e.g.,  McEwan  and  Phillips,  1975) 

b.  Satellite  Observations 

Satellite  airglow  observations  offer  an 
advantage  over  ground  observations  in  that  there  is  almost 
no  scattered  sunlight  to  obscure  the  actual  airglow 
emissions.  An  additional  advantage  lies  in  the  geometry  of 
satellite  limb  viewing  (see  e.g.,  Liou,  1980).  This  makes 
possible  the  calculation  of  the  column  density  of  an 
emitting  species.  Through  an  "onion-skin"  inversion  method 
the  number  density  of  the  emitting  species  can  then  be 
calculated  (see  e.g.,  Cleary,  1985;  Hickman  et  al.,  1977). 

2.  Radio  Sounding  of  the  Ionosphere 

Electron  density  profiles  of  the  ionosphere  can  be 
measured  using  a  variable  frequency  radar  which  measures  the 
time  delay  between  the  transmission  of  an  RF  pulse  and  its 
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echo  from  a  reflecting  level  within  the  ionosphere  (see 
e.g.,  Bauer, 1973).  Such  a  device  is  called  an  ionosonde . 

The  requirement  for  reflection  of  an  electromagnetic 
wave  from  an  ionized  medium  is  that  the  magnitude  of  the 
complex  propagation  vector,  K ,  be  imaginary.  This  require¬ 
ment  is  met  when 

[e-~\  >  1.  (2.1) 

e0  moo 


The  value  of  to  for  which  equation  (2.1)  is  identically  one 
is  <Op,  the  plasma  frequency.  The  index  of  refraction,  n,  in 
the  ionized  medium  is  given  by 


n 


a 


to 


1  - 


80 . 6 [ e ] 


(2.2) 


The  requirement  for  reflection  of  an  electromagnetic  wave  at 
normal  incidence  is  n  =  0.  Thus  the  maximum  frequency  for 
normal  reflection  is  given  by 


f2  -  80 . 6 [ e ] 


(2.3) 


where  f  is  measured  in  Hz  and  [e]  is  the  electron  density  in 
electrons/m2.  Thus  by  varying  the  frequency  it  is  possible 
to  map  out  the  electron  density  profile.  However, 
examination  of  Figure  2.1  shows  that  the  electron  density 
above  about  250  km  cannot  be  measured  by  ground-based 
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ionosondes  because  the  frequencies  needed  to  probe  this 
region  will  already  have  been  reflected  at  a  lower  altitude. 

In  order  to  obtain  electron  density  profiles  for  the 
region  above  250  km,  ionosondes  have  been  placed  on 
satellites.  Top-side  sounding  is  a  very  valuable  technique 
for  obtaining  electron  density  profiles  in  the  upper 
F2-region.  Unfortunately,  most  of  the  measurements  taken 
from  top-side  ionosondes  do  not  provide  full  profiles  down 
to  the  peak  altitude  so  that  it  is  difficult  to  merge  this 
data  with  that  of  ground-based  ionosondes. 
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III.  CHEMICAL  REACTIONS  IN  THE  THERMOSPHERE 


A  great  variety  of  atoms  and  molecules  exist  in  the 
thermosphere.  As  will  be  discussed  later,  much  of  this 
molecular  soup  is  subject  to  reaction  with  solar  radiation. 
The  resulting  photolytically  generated  species  undergo  many 
different  chemical  processes  with  the  surrounding 
atmosphere.  This  chapter  will  examine  these  processes  and 
introduce  the  method  used  for  analyzing  the  rapidly  changing 
composition  and  structure  of  the  thermosphere. 

A.  PRINCIPLES  OF  PHOTOCHEMISTRY 

Photochemistry  is  concerned  with  chemical  reactions  that 
are  initiated  by  the  absorption  of  solar  radiation.  The 
absorption  process  is  written  as 

A  +  hv  =►  A* ,  (3.1) 

* 

where  A  is  an  excited  energy  state  of  molecule  A.  The 
excited  state  could  be  a  rotational  or  vibrational  level,  an 
electronic  state,  or  it  could  represent  the  decomposition  of 
the  molecule  into  two  or  more  fragments.  This  latter 
process  is  called  photodissociation  or  photolysis  if  the 
products  are  all  neutral  atoms  or  molecules,  and  is  called 
photo  ionization  if  the  products  include  an  ion  and  one  or 
more  electrons. 
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Once  excited  molecules  or  atoms  are  produced,  they 
undergo  secondary  chemical  reactions,  forming  species  that 
would  not  exist  in  a  normal  state  of  local  ^rmal 

equilibrium  ( LTE ) .  Therefore  these  photochemical  processes 
are  non-LTE  mechanisms  and  must  be  quantitatively  analyzed 
by  studying  the  rate  at  which  reaction  occurs  (see  e.g., 
Chamberlain, 1978 ) . 

1.  Photodissociation 

The  process  by  which  the  absorption  of  a  photon 
leads  to  disruption  of  the  absorbing  molecule  is  called 
photodissociation  or  photolysis.  The  atoms  produced  in  this 
manner  are  highly  reactive.  Photolysis  processes  which  take 
place  in  the  photochemical  model  used  for  this  thesis  are: 

NO  +  hv  *  O  +  N(4S) ,  (3.2) 

N2  +  hv  *  N(4S)  +  N(4S),  (3.3) 

and 

N2  +  hv  N(2D)  +  N(2D).  (3.4) 

Notice  that  the  last  two  equations  involve  the 
photodissociation  of  molecular  nitrogen  into  two  different 
electronic  states  of  elemental  nitrogen.  The  extent  to 
which  these  two  reactions  occur  is  given  by  their  branching 
ratios .  For  the  above  two  reactions,  the  branching  ratio 
for  each  has  been  determined  to  be  0.5,  indicating  that  half 
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of  the  reactions  which  occur  will  produce  N(4S)  and  that  the 

2 

other  half  will  result  in  N(  D)  production  (see  e.g., 
Heiklen,  1976 )  . 

a.  Photodissociation  Rates 

There  are  three  factors  influencing  the  rate  of 
photodissociation  :  molecular  concentration,  photon  flux, 
and  quantum  yield.  The  most  obvious  factor  must  be  the 
amount,  or  concentration,  of  molecule  A.  Concentrations  of 
a  substance  per  unit  volume  are  denoted  by  enclosing  the 
substance  in  square  brackets.  [A]  is  the  concentration  of 
substance  A  per  unit  volume. 

Another  factor  influencing  the  rate  is  the 
number  of  photons  available  to  cause  photodissociation.  The 
flux  of  solar  radiation  outside  the  atmosphere  is  simply  a 
function  of  the  brightness  of  the  Sun  and  the  distance  from 
the  Sun.  However,  as  the  radiation  penetrates  the  Earth's 
atmosphere  the  flux  of  photons  is  attenuated  due  to 
absorption  by  atmospheric  gasses.  At  high  altitudes  the 
atmospheric  density  is  low  and  therefore  little  absorption 
takes  place.  As  the  altitude  decreases  the  density 

increases,  and  absorption  becomes  important. 

The  efficiency,  or  quantum  yield  of  the  reaction 

* 

given  by  equation  (3.1)  is  defined  as  the  number  of  A 
molecules  produced  for  each  absorbed  photon.  Quantum  yield 
varies  with  wavelength  and  is  zero  longward  of  the 
ionization  wavelength. 
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A  common  way  of  expressing  the  photon  flux  at  a 
particular  altitude  z  is  in  terms  of  the  exoatmospher ic  flux 
and  the  optical  depth.  Optical  depth,  x,  is  defined  as 

t  -  1 1  ]dz '  ,  (3.5) 

i  J’ 

I.  1* 

where  a  is  the  absorption  cross  section  for  the  ic 
constituent  at  a  specific  wavelength,  and  [i]  is  the 
concentration  of  that  constituent.  The  photon  flux  is  then 
given  by 

*(z)  -  f_e-T'sec(x),  (3.6) 

where  *(z)  is  the  photon  flux  at  altitude  z,  is  the  exo- 
atmospheric  photon  flux,  x  is  optical  depth,  and  x  is  the 
solar  zenith  angle.  Note  that  x*sec(x)  is  the  optical 
pathlength  through  the  atmosphere. 

2.  Photoionization 

When  the  products  of  equation  (3.1)  include  molec¬ 
ular  or  atomic  ions  and  one  or  more  electrons,  the  process 
is  called  photoionization.  Schematically  this  is  written 

A  +  hv  4  a+  +  e*.  (3.7) 

In  this  process,  the  energy  carried  in  by  the  photon  is 
converted  into  ionization  energy  plus  kinetic  energy  of  the 
electron.  These  energetic  electrons  are  called  photo- 
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* 

electrons,  (e  ),  and  enter  into  several  different  energy 
transfer  reactions  with  neutral  gasses. 

The  major  source  of  ionizing  radiation  in  the 
thermosphere  is  solar  radiation.  The  only  wavelengths  which 
penetrate  to  the  E-layer  are  extreme  ultraviolet  (EUV), 
X-rays,  and  Lyman-a. 

a.  Photoionization  Rates 

The  rate  at  which  photoionization  takes  place 
for  each  constituent  is  described  by  its  ionization  rate. 
For  a  given  constituent,  the  ionization  rate  is  given  by 

qi(z)  =  J  <VX)‘*z(X)  dX'  (3.8) 

where  a  is  the  ionization  cross  section  and  $  is  the  solar 

z 

flux  at  altitude  z.  (Torr  et  al.,  1979) 

B.  DIFFUSION  AND  THE  CONTINUITY  EQUATION 

It  is  possible  to  construct  an  equation  showing  the  time 
variation  of  a  particular  constituent's  concentration  once 
the  rate  coefficients  for  the  production  and  loss  mechanisms 
are  known.  This  equation  is  called  the  continuity  equation 
and  contains  terms  relating  the  chemical  creation  and 
destruction  of  the  constituent,  as  well  as  terms  describing 
the  flux  due  to  physical  transport  processes.  There  are 
four  types  of  physical  transport  mechanisms:  molecular 
diffusion,  ambipolar  diffusion,  eddy  diffusion,  and  thermal 


1.  Molecular  and  Ambipolar  Diffusion 

When  the  concentration  gradient  of  a  constituent  in 
the  atmosphere  differs  from  its  equilibrium  value,  the 
resulting  motion  is  called  molecular  diffusion.  For  an 
isothermal  atmosphere,  this  can  be  written 


«. 

1 


3[i] 

dz 


(3.9) 


where  is  the  flux,  is  the  molecular  diffusion  coef- 

u 

ficient,  and  is  the  scale  height  for  the  ic  constituent. 

Diffusion  which  results  from  a  concentration 
gradient  of  a  charged  constituent  in  a  partially  ionized 
gas,  or  plasma,  in  which  the  constituent  is  influenced  by 
the  electrostatic  forces  of  the  plasma  as  a  whole  is  called 
ambipolar  diffusion.  Although  electrons  in  the  thermosphere 


have  a 

greater  diffusion 

coefficient 

than 

the  ions, 

electrostatic  attraction  between  electrons 

and  positive  ions 

prevents 

any  large  scale 

separation . 

Thus 

the  plasma 

diffuses 

as  a  whole,  with 

ambipolar  diffusion 

coefficient 

D  .  (see  e.g.,  McEwan  and  Phillips,  1975) 

2.  Eddy  Diffusion 

Eddy  diffusion  is  the  result  of  turbulent  mixing  in 
the  presence  of  a  composition  gradient.  Whereas  molecular 
diffusion  results  from  the  average  velocity  of  the 
constituent  being  different  from  that  of  others,  eddy 
diffusion  occurs  whether  or  not  there  is  any  difference  in 
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average  molecular  velocity  between  the  constituents.  Eddy 
diffusion  can  be  expressed  similarly  to  molecular  diffusion 
by  an  equation  analogous  to  (3.8).  The  coefficient  of 
molecular  diffusion  is  replaced  by  the  eddy  diffusion 
coefficient  K. 

3.  Thermal  Diffusion 

Thermal  diffusion  is  due  to  large  temperature 
gradients  in  a  gas  mixture.  It  can  be  properly  accounted 
for  only  in  a  rather  complex  manner.  However,  since  the 
thermal  diffusion  factor  ai  is  proportional  to  the 
difference  in  molecular  weights  of  the  constituent  gasses, 
thermal  diffusion  is  important  in  the  atmosphere  only  for 
light  gasses  such  as  hydrogen  and  helium.  Because  the  model 
used  in  this  thesis  is  primarily  concerned  with  the 
chemistry  of  nitric  oxide,  thermal  diffusion  can  be 
neglected . 

4.  The  Continuity  Equation 

Physical  transport  and  chemical  processes  act  on 
each  constituent  of  the  atmosphere  to  produce  a  variation  in 
concentration  with  time.  The  equation  which  governs  these 
processes  is  the  mass  continuity  equation 


Pi 


L.[i] . 


(3.10) 


16 


Again,  [i]  is  the  concentration  at  time  t  and 
altitude  z,  and  L^[i]  are  chemical  production  and  loss 
rates,  and  is  the  mean  vertical  velocity  due  to  transport 
processes.  Since  [i]*V^  is  the  vertical  flux  #^,  the  second 
term  in  equation  (3.9)  may  be  written 


3  ♦ 

‘  TT-  l3-U) 

For  a  neutral  constituent  with  low  concentration,  may  be 
written  as  the  sum  of  molecular  and  eddy  diffusion.  Thus 


*. 

x 


'Di 


'  3  l  i  1  .  t 


9  z 


jj 

T 


JJ  4.  JJJ)  _  K  ra  [  i  )  +  [  i  ]  3T  .  I  i  n 

9z  +  H  J  *  rtr  +  T  ~9z  +  “H~  y 

l  9  V 


(3.12) 


where  the  middle  terms  are  due  to  the  dependence  of  scale 
height  on  temperature  gradient,  and  K  are  the  molecular 
and  eddy  diffusion  coefficients  respectively,  T  is  the 
absolute  temperature,  Hi  the  scale  height  for  the  ith 
constituent,  and  H  the  scale  height  for  a  constituent 
having  the  mean  molecular  mass.  (Colgrove  et  al.,  1966) 
a.  Time  Constants 

The  continuity  equation  includes  chemical  and 
dynamical  processes.  The  dominance  of  one  process  over 
another  can  be  determined  by  examining  the  relative 
magnitude  of  the  characteristic  times,  or  time  constants,  of 
the  competing  processes  (Goody  and  Walker,  1972).  If  two 


processes  are  competing  and  one  acts  more  rapidly  than  the 


other,  the  faster  one  will  be  more  effective.  The  dominant 
process  will  be  the  one  with  the  shorter  time  constant  (see 
e.g.,  Wallace  and  Hobbs,  1977). 

By  performing  dimensional  analysis  on  the 
continuity  equation,  the  time  constants  for  chemistry,  x  , 
eddy  diffusion,  xR,  and  molecular  diffusion,  xD,  can  be 
approximated  (Cleary,  1985).  The  distribution  of  [i]  is 
determined  principally  by  chemistry  if  xc  <  rR  and  xc  <  td, 
and  by  dynamics  if  xR  <  xc  and  xD  <  tc. 

A  comparison  of  the  lifetimes  for  molecular 
diffusion,  eddy  diffusion,  and  chemical  destruction  was  done 
by  Cleary  (1985)  for  the  different  species  used  in  that 
model.  The  comparison  showed  that  dynamical  processes  are 

4 

only  important  for  two  species,  NO  and  N(  S).  For  the 
remaining  species  (N(2D),  N0+ ,  02+,  N-,+  ,  0+  •  and  0(2D)  +  ), 
chemical  processes  are  sufficiently  fast  that  transport  may 
be  neglected  (Cleary,  1985).  Ambipolar  diffusion  for  0+  was 
later  added  to  the  model. 

5.  Solutions  to  the  Continuity  Equation 

For  constituents  which  have  sufficiently  short 
chemical  time  constants,  the  dynamical  processes  of 
diffusion  may  be  neglected.  Thus  the  flux  term  in  equation 
(3.10)  vanishes  and  the  continuity  equation  reduces  to 

-  P  -  L[i]  (3.10a) 
for  N(2D),  NO+,  02\  N2+,  0+,  and  0(2D)+. 
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For  those  constituents  which  have  short  dynamical 

4 

time  constants  (i.e.,  NO  and  N(  S)),  the  flux  term  in  the 
continuity  equation  must  be  retained.  When  the  expression 
for  *  given  by  equation  (3.12)  is  substituted  into  equation 
(3.10)  and  differentiated,  a  parabolic  differential  equation 
results.  This  equation  can  be  written 


3[i) 

3t 


(3.13) 


The  coefficients  A^ ,  ,  C^,  and  Ei  in  equation 
(3.13)  are  given  in  Table  I.  Equations  (3.10a)  and  (3.13) 
cannot  be  solved  analytically  and  must  be  solved  by 
numerical  methods.  In  the  model  used  for  this  thesis,  the 
method  is  the  Crank-Nicolson  finite  differencing  technique. 
(Cleary,  1985  App.  A) 


C.  CHEMICAL  REACTIONS  IN  THE  LOWER  THERMOSPHERE 

Many  types  of  chemical  reactions  are  continually  taking 
place  in  the  atmosphere.  Products  of  these  reactions  include 
photolytically  generated  atoms,  radicals  excited  molecules 
and  ions.  As  was  seen  in  previous  sections,  photoionization 
and  photodissociation  are  the  primary  photochemical 
reactions.  Reactions  involving  the  products  of  these  two 
types  of  reactions  with  themselves  and  with  the  stable 
atmospheric  gasses  determine  in  great  part  the 
characteristics  of  the  atmosphere. 
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TABLE  I 

COEFFICIENTS  TO  THE  PARABOLIC  DIFFERENTIAL  EQUATION 


The  appendix  summaries  the  types  of  chemical  processes 
which  are  important  in  the  execution  of  the  one  dimensional 
model  used  in  this  thesis  for  the  calculation  of  nitric 
oxide  in  the  thermosphere.  Of  course,  many  other  reactions 
also  occur,  but  it  is  the  chemistry  of  nitric  oxide  which  is 
germane  to  this  investigation. 
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1.  Nitric  Oxide  Chemistry 

Direct  sources  of  nitric  oxide  in  the  E-region  are 

4 

the  reactions  of  atomic  nitrogen  (ground  state  N(  S)  and 

2  2 

excited  states  N(  D)  and  N(  P))  with  molecular  oxygen. 
These  three  reactions  (R28,  R29,  and  R39  of  Appendix)  do  not 
equally  contribute  to  the  production  of  NO.  In  steady 

state,  the  NO  production  rate  of  each  reaction  is  found  by 
multiplying  the  rate  coefficient  by  the  product  of  the 
concentrations  of  the  reactants.  Figure  3.1  shows  some 

typical  daytime  nitrogen  density  profiles.  Figure  3.2  shows 
the  corresponding  NO  production  rates  from  the  three 


different 

atomic 

nitrogen  states,  as  well 

as 

the 

rate 

coefficient 

(k) 

for  each  reaction.  It  can 

be 

seen 

from 

Figure  3.2 

that 

the  principal  source  of  NO 

at 

110 

km  is 

2 

N (  D)  via  the  reaction 

N(2D)  +  02  *►  NO  +  0  (3.14) 

4 

The  small  contribution  from  N(  S)  is  due  to  the  temperature 

dependence  of  the  rate  coefficient,  and  the  small 

2 

contribution  from  N(  P)  is  because  of  its  low  concentration. 

2 

Sources  of  N(  D)  are  from  the  dissociative 
recombination  of  N2+  and  N0+  (R19  and  R27),  the  dissociative 
excitation  by  photoelectrons  of  N2  (R12),  the  photo¬ 
dissociation  of  N2  (R15),  and  the  ion-atom  interchange 
between  and  0  (R17).  In  a  method  analogous  to  that  used 

for  NO  production,  it  can  be  shown  that  the  two  principal 
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Altitude  (km) 


2 

sources  of  N(  D)  are  via  R27  and  R12.  These  equations  are 
stated  here  for  ease  of  reference 

N0+  +  e  •»  0  +  N(2D)  R27 

N2+  +  e*  *  N(2D)  +  N(2D)  +  e  R12 

The  two  loss  processes  for  nitric  oxide  in  the  E-region  are 

4 

via  the  bimolecular  reaction  with  N(  S),  (R33),  and  the 


charge  exchange  with  02+ ,  (R22) 

N(4S)  +  NO  4  N2  +0 

R33 

02+  +  NO  ■»  N0+  +  02 

R22 

22 


Altitude  (km) 


Figure  3.2.  NO  Production  Rates 


2 

Thus,  while  N(  D)  is  a  source  of  nitric  oxide  at  100  km, 

4 

N (  S)  is  both  a  source  and  sink.  An  analysis  of  these 
production  and  loss  processes  as  they  pertain  to  the 
photochemical  model  is  presented  in  Chapter  V. 
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IV.  THE  PHOTOCHEMICAL  MODEL 


A.  PHOTOCHEMICAL  MODELING 

Photochemical  modeling  is  a  useful  tool  for  gaining  an 
understanding  of  the  minor  constituents  in  the  thermosphere. 
By  comparing  a  theoretical  model  with  observations  of  some 
constituents  for  which  data  are  available,  one  can  determine 
the  validity  of  the  model  and  better  understand  physical  and 
chemical  mechanisms  occurring  in  the  thermosphere. 

"Photochemical  models  are  sensitive  to  various  input 
parameters  such  as  solar  fluxes,  photoabsorption 
cross-sections,  chemical  reaction  schemes  and  reaction  rate 
constants,  and  transport  parameters  (e.g.,  eddy  diffusion 
coefficient)  as  well  as  to  the  method  of  numerical  analysis 
and  boundary  conditions."  (Shimazaki,  1985).  Additionally, 
models  may  be  time  dependent  or  time  independent.  Time 
independent  models  are  usually  used  to  arrive  at  a 
steady-state  condition.  A  steady-state  model  does  not 
represent  the  actual  state  for  the  constituents  whose 
concentrations  change  as  a  function  of  time.  Causes  for 
these  changes  would  be  solar  flux,  geomagnetic  Ap  index,  and 
any  temperature  dependent  reaction  rate  since  the  temp¬ 
erature  is  also  time  dependent.  Time  dependent  models 
attempt  to  account  for  these  temporal  variations. 
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B.  THE  TINE  DEPENDENT  PHOTOCHEMICAL  MODEL 
1.  Numerical  Technique 

In  Chapter  III  it  was  shown  that  the  three  phenomena 

which  must  be  understood  in  order  to  model  an  atmospheric 

constituent  properly  are  chemistry,  eddy  diffusion,  and 

molecular  diffusion.  Additionally,  it  was  shown  that  the 

continuity  equation  reduces  to  equation  (3.10a)  for  species 

which  do  not  require  transport  processes  to  be  considered. 
4 

For  NO  and  N(  S)  transport  processes  must  be  included  and 
the  continuity  equation  becomes  a  parabolic  differential 
equation  (equation  (3.13)). 

Equation  (3.10a)  can  be  solved  by  the  following 
semi-implicit  finite  difference  equation 


At  *  *t 

LtU 

]t+l 

9 

(4.1) 

which  is  solved  for  [i]^+^,  where  [i 

]t  +  l 

is 

the 

density  of 

the  i*"*1  constituent  at  time  interval 

t+1 

and 

is 

a  function 

of  [ i ] t . 

Equation  (3.13)  is  first 

order 

in 

the  time 

derivative  and  second  order  in  the  spatial  derivative  and 
therefore  requires  three  boundary  conditions:  one  in  time 
and  two  in  space.  Cleary  (1985)  has  shown  and  Siskind 
(1988)  has  verified  that  the  spatial  boundary  conditions  are 
satisfied  by  specifying  conditions  at  the  upper  and  lower 
altitude  limits  of  the  model.  Several  different  spatial 
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conditions  can  be  used  to  implement  these  boundary 

conditions.  The  ones  used  by  the  model  require  that  there 

4 

be  zero  flux  of  both  NO  and  N(  S)  at  the  upper  boundary, 
and  at  the  lower  boundary  the  condition  for  NO  is  that  there 
is  a  local  minimum  in  the  density  profile,  while  for  N(^S) 
the  requirement  is  that  it  be  in  a  state  of  photochemical 
equilibrium. 

When  equation  (3.13)  is  cast  into  its  finite 
difference  form  using  the  Crank-Ni colson  method,  a  series  of 
coupled  algebraic  equations  are  produced  which  can  be 
expressed  as  a  tri-diagonal  matrix  equation.  The  density 
profile  for  each  time  increment  is  obtained  using  the  Thomas 
algorithm  of  matrix  inversion.  (Von  Rosenberg,  1969;  Cleary, 
1985) 

2.  Inputs  to  the  Model 

The  Cleary  time  dependent  photochemical  model  is  run 
in  three  distinct  steps.  The  first  step  uses  the  MSIS-83 
program  (Hedin,  1983)  to  calculate  the  neutral  atmosphere. 
The  second  step  uses  a  modification  of  the  NRL  photoelectron 
code  PEGFAC  (NRL  Memorandum  Report  5004,  1982)  to  calculate 
photoelectron  impact  rates.  Siskind  (1988)  modified  the 
PEGFAC  code  to  extend  the  energy  grid  so  as  to  include  soft 
X-rays  up  to  666  eV  (18.62  A).  The  existing  energy  grid  was 
also  updated  to  reflect  more  recent  data.  The  third  step  is 
the  actual  one  dimensional  model.  The  model  calculates  the 
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time  development  of  each  species  which  enters  into  the 
chemistry  of  nitric  oxide. 

a.  Neutral  Atmosphere 

The  first  step  of  the  photochemical  model  uses 
the  MSIS-83  (Hedin,  1983)  code  to  calculate  the  neutral 
atmosphere  for  each  model  "hour"  of  every  day  for  which  the 
model  is  run.  Results  are  stored  as  183  element  column 
vectors.  Each  "hour"  five  of  these  vectors  are  generated. 
The  first  contains  information  on  the  geophysical  parameters 
used.  These  are  (1)  the  exoatmospheric  temperature,  (2) 
solar  zenith  angle,  (3)  the  F10.7  value  for  the  previous 
day,  (4)  the  90  day  average  F10.7  value,  (5)  geomagnetic  Ap 
index,  (6)  latitude,  (7)  longitude,  and  (8)  the  X-ray 
scaling  factor,  FRACXR .  The  second  vector  contains  the 
temperature  profile  for  altitudes  from  86  to  450  km  in  2  km 
increments  (hence  183  elements).  The  third  fourth  and  fifth 
column  vectors  contain  the  density  profiles  of  O,  C^,  and 
N2,  respectively.  Thus  each  "hour"  is  a  183  x  5  array. 
Each  day  produces  24  such  arrays. 

b.  Photoelectron  Impact  Rates 

The  calculation  of  photoelectron  impact  rates 
was  done  using  the  computer  program  PEGFAC  ( NRL  Memorandum 
Report  5004,  1982).  Recently,  Siskind  (1988)  revised  the 

PEGFAC  code  to  account  more  properly  for  the  higher  energy 
photoelectrons  produced  by  soft  X-rays.  He  expanded  the 
original  37  EUV  intervals  from  50-1050  A  (Torr  et  al . ,  1979) 
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to  include  an  additional  20  intervals  which  cover  the  soft 
X-ray  flux  from  18-50  A.  The  fluxes  used  were  measured  in 
one  of  Hinteregger '  s  rockets  and  are  available  from  the 
National  Space  Science  Center  in  a  file  labelled  SC#2lREFW 
( Siskind,  1988 )  . 

EUV  fluxes  are  calculated  by  performing  a  linear 

interpolation  between  the  appropriate  solar  flux  periods 

reported  by  Torr  et  al.  {1979)  based  on  the  value  of  F10.7 
-22  -2  -1 

(xlO  W  m  Hz  ).  Thus  the  EUV  flux  is  scaled  according 
to  F10.7.  Since  only  one  measurement  of  X-ray  flux  is 
available,  interpolation  is  not  possible.  However,  in  order 
to  explain  the  observed  solar  cycle  variations  of  NO,  Barth 
et  al.  (1988)  suggested  that  the  soft  X-ray  flux  exhibits  a 
larger  variation  with  solar  activity  than  EUV.  For  this 
reason,  Siskind  added  a  multiplicative  X-ray  flux  scaling 
factor,  FRACXR,  to  the  PEGFAC  code.  The  X-ray  fluxes  used 
by  the  model  are  scaled  by  the  multiplicative  factor  FRACXR. 
It  is  the  goal  of  this  thesis  to  determine  the  relationship 
between  FRACXR  and  F10.7. 

PEGFAC  uses  the  flux  data  to  calculate  the 
photoionization  rates  of  N£,  O2,  and  0,  and  the  photo¬ 
dissociation  rates  of  Nj  for  each  model  "hour"  of  every  day 
for  which  the  model  is  run.  As  with  MSIS-83,  PEGFAC  stores 
the  results  as  183  element  column  vectors.  For  each  "hour", 
four  of  these  vectors  are  generated.  The  first  three 
contain  the  photoionization  rates  for  N2 ,  Oj,  and  0, 


respectively.  The  fourth  contains  the  Nj  photodissociation 
rates.  Thus  each  hour  is  a  183  x  4  array.  Each  day 
produces  24  such  arrays. 

c.  Direct  Inputs  to  the  Model 

The  neutral  atmosphere,  photoionization  rates 
for  N2,  O2 ,  and  0,  and  the  photodissociation  rates  for  N2 
are  calculated  exterior  to  the  model.  The  data  files 
created  by  these  two  procedures  are  used  as  inputs  to  the 
model  program.  The  density  profiles  for  the  minor  species 
are  set  to  a  generic  initial  condition,  and  the  model  is 
allowed  to  run  in  a  time  independent  mode  for  a  sufficient 
length  of  time  to  allow  this  initial  minor  atmosphere  to 
reach  the  steady  state  conditions  which  correspond  to  the 
starting  date  of  the  run.  The  model  "solves"  the  continuity 
equation  for  each  of  the  species  involved  in  the  nitric 
oxide  chemistry  every  15  minutes  using  the  numerical  methods 
previously  described.  When  the  sun  is  within  15  degrees  of 
the  horizon  the  atmosphere  is  rapidly  changing,  therefore 
the  time  increment  of  the  model  switches  to  five  minutes. 
As  with  MSIS-83  and  PEGFAC,  the  model  stores  the  density 
profile  for  each  minor  constituent  in  a  183  element  column 
vector . 

3.  Outputs  From  the  Model 

The  program  is  written  to  calculate  the  density 
profile  for  the  minor  species  at  15  minute  intervals,  or  at 
five  minute  intervals  if  the  sun  is  within  15  degrees  of  the 
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horizon.  For  the  purpose  of  this  thesis,  the  output  has 
been  structured  to  extract  the  nitric  oxide  density  at  110 
km  for  each  day  at  1500  local  time.  This  is  because  the  SME 
satellite  data  was  taken  at  an  altitude  of  110  km  at  1500 
local  time  each  day.  The  output  is  stored  in  a  vector 
containing  N  elements,  where  N  is  the  number  of  days  for 
which  the  model  was  run. 
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V.  ANALYSIS 


Chapters  III  and  IV  introduced  the  theory  and  modeling 
technique  that  is  used  in  this  thesis.  This  chapter  will 
present  the  experimental  approach  used  to  determine  the 
relationship  between  FRACXR  and  F10.7. 

A.  SHE  OBSERVATIONS 

Data  from  the  Solar  Mesosphere  Explorer  satellite  was 
provided  by  Dr.  Charles  A.  Barth  of  the  Laboratory  for 
Atmospheric  and  Space  Physics,  University  of  Colorado 
(telephone  conversation  17  November,  1988).  Equatorial 
latitude  nitric  oxide  data  was  used  for  this  thesis.  This 
data  is  the  result  of  averaging  measurements  from  two 
satellite  orbits  between  20°  north  latitude  and  20°  south 
latitude.  The  first  orbit  took  measurements  from  20°-60° 
east  longitude,  while  the  second  were  between  105°-147°  west 
longitude.  All  measurements  were  taken  at  an  altitude  of 
110  km  at  1500  local  time.  The  110  km  nitric  oxide  density 
for  the  period  January  6,  1982  to  April  19,  1985  is  shown  in 
Figure  5.1. 

1.  Analysis  of  Original  Model 

Figure  5.2  shows  a  comparison  between  SME  data  and 
model  calculations  prior  to  the  addition  of  Siskind's  FRACXR 
parameter.  It  clearly  shows  that  the  model  was  inadequate 
in  predicting  both  the  absolute  magnitude  and  variations  in 
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1982  1983  1984  1985 

Figure  5.1.  Nitric  oxide  density  at  110  km. 

SME  measurements  from  January  6,1982  to 
April  19,  1985. 

the  nitric  oxide  density  which  were  observed  by  the  SME 

satellite.  The  reason  for  this  inadequacy  can  be  explained 

by  taking  a  closer  look  at  the  relevant  photochemistry. 

Originally,  an  increase  in  the  EUV  flux  would  increase  the 

2 

rate  at  which  the  production  of  N(  D)  occurs  (see  R12  and 

R27  in  Appendix).  However,  this  increase  in  EUV  flux  would 

also  cause  one  of  the  principal  loss  reactions  for  nitric 

oxide  (R22)  to  increase  its  rate  due  to  an  increase  in  the 

2 

ionized  molecular  oxygen  density.  The  increase  in  N(  D) 
(and  hence  NO)  was  thereby  partially  offset  by  the  increased 
loss  of  NO  via  R22.  A  process  was  sought  which  would  allow 
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Data 

I 

Figure  5.2.  Comparison  between  SME  observations 
and  the  photochemical  model  prior  to  FRACXR. 

2 

the  model  to  increase  the  production  of  N(  D)  without 

causing  a  corresponding  decrease  in  NO  density  via  the 

loss  reaction  R22.  This  led  to  the  modification  of  the 

PEGFAC  code  by  Siskind  as  discussed  in  Chapter  IV. 

The  reason  behind  the  incorporation  of  the 

multiplicative  X-ray  scaling  factor  is  now  clear.  Soft 

X-rays  which  penetrate  to  the  E-region  photoionize  the  major 

constituents  and  produce  energetic  photoelectrons  as 

discussed  in  Chapter  III.  The  increase  in  photoelectrons 

2 

drives  the  production  of  N(  D)  via  R12  and  R17,  but  does  not 
affect  R22.  In  order  for  this  mechanism  to  work,  the 
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variation  in  solar  X-radiation  must  be  greater  than  the 
variation  in  solar  EUV  radiation.  (Barth  et  al.,  1988) 

B.  DETERMINATION  OF  FRACXR 

initial  estimates  of  the  magnitude  of  FRACXR  were 
provided  by  Siskind  (telephone  conversation  November  9, 
1988).  He  suggested  that  FRACXR  take  on  values  near  one  for 
conditions  of  solar  minimum  and  50  for  solar  maximum.  This 
provided  a  starting  point  for  the  determination  of  the 
relationship  between  FRACXR  and  solar  cycle. 

1.  Assumptions 

Two  major  assumptions  were  made  for  the  purposes  of 
simplifying  the  data  analysis  and  reducing  computer  time. 
The  first  assumption  is  that  auroral  electrons  do  not 
influence  the  nitric  oxide  chemistry.  This  assumption  is 
based  on  the  fact  that  all  computations  and  observations 
were  at  equatorial  latitudes.  It  is  justified  by  the  work 
of  Siskind  (1988),  who  showed  that  for  equatorial  latitudes 
there  is  little  enhancement  of  NO  density  due  to  an  auroral 
storm. 

The  second  assumption  was  that  the  computation  of 
nitric  oxide  density  would  be  independent  of  the  longitude 
at  which  the  model  was  run.  This  allowed  the  model  to  be 
run  at  only  one  longitude  for  the  computation  of  each  NO 
density  time  sequence,  instead  of  averaging  two  sets  of 
measurements  as  was  done  for  the  SME  data.  To  justify  this 
assumption,  the  model  was  run  at  both  40  degrees  east 


longitude  and  126  degrees  west  longitude  for  the  same  five 
day  period  in  1982.  Figure  5.3  shows  the  result  of  this 
procedure.  The  two  cases  deviated  from  each  other  by  less 
than  0.2%.  Based  on  this  result,  all  subsequent  model  runs 
were  done  only  at  40  degrees  east  longitude. 


Date 

Figure  5.3.  Model  (NO]  calculations  for 
two  different  longitudes. 


2.  Experimental  Approach 

The  SME  data  was  searched  to  locate  periods  where 
the  value  of  F10.7  remained  constant  for  at  least  three 
days.  Four  such  periods  were  initially  found  corresponding 
to  F10.7  (xlO-22  Wm-2  Hz-1)  values  of  70,  90,  128,  and  201. 
The  model  was  then  run  with  various  values  of  FRACXR  for 


35 


each  of  the  four  periods.  Figures  5.4  through  5.7  show  the 
results  of  these  runs.  This  initial  data  is  summarized  in 
Table  II  and  shown  graphically  in  Figure  5.8.  Uncertainties 
are  estimated  from  the  data  directly  without  measuring  them 
explicitly  since  only  one  set  of  calculations  is  possible 
each  time  the  model  is  run  (see  e.g.f  Bevington , 1969 ) . 


SUMMARY 

TABLE  II 

OF  INITIAL  FRACXR  DATA 

DATE 

F10.7 

START 

END 

(xl0~22  Wm-2Hz_1) 

FRACXR 

84300 

84303 

70 

3  ±  2 

83297 

83300 

90 

6  ±  1 

83185 

83188 

128 

22  ±  2 

82169 

82171 

201 

24  ±  1 

It  can  be  seen  from  Figure  5.8  that  the  first  three 
data  points  are  nearly  linear,  while  the  point  corresponding 
to  an  F10.7  value  of  201  lies  significantly  lower  than  would 
be  expected  using  a  linear  relationship  based  on  the  first 
three  points.  To  resolve  this  uncertainty,  a  fifth  data 
point  was  taken  for  an  F10.7  value  of  242.  Although  this 
value  was  constant  for  only  two  days,  the  FRACXR  value  at 
this  point  would  determine  if  a  linear  relationship  was 
probable.  The  result  of  this  run  is  shown  in  Figure  5.9. 
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Figure  5.4.  Model  INOl  densities  for  F10.7  -  70. 
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iqure  5.5.  Model  (NO]  densities  for  F10.7 
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3.  Fitting  the  Data 


The  addition  of  the  fifth  data  point  at  F10.7  -  242 
is  reflected  in  Table  III,  a  summary  of  the  final  FRACXR 
data.  The  problem  now  is  to  find  a  function  (or  combination 
of  functions)  which  best  fits  the  data.  This  function  can 
then  be  used  to  determine  the  value  of  FRACXR  on  a  daily 
basis  for  input  to  the  model. 


SUMMARY 

TABLE  III 

OF  FINAL  FRACXR  DATA 

DATE 

START  END 

F10.7 

(xlO-22  Wm-2Hz-1 ) 

FRACXR 

84300 

84303 

70 

3  ±  2 

83297 

83300 

90 

6  ±  1 

83185 

83188 

128 

22  ±  2 

82169 

82171 

201 

24  ±  1 

82342 

82343 

242 

27  ±  1 

Several  fitting  functions  were  tested.  The  criteria 

for  determining  the  best  fit  was  that  function  which 
2 

minimized  X  and  was  monotonically  increasing  in  FRACXR. 

2 

"The  definition  of  X  as  the  ratio  of  the  estimated  variance 
2  2 

s  to  the  parent  variance  a  ...makes  it  a  convenient  measure 
of  the  goodness  of  fit."  (Bevington,  1969).  A  weighted 
least-squares  fit  to  a  third  degree  polynomial  is  shown  in 
Figure  5.10,  and  a  weighted  least-squares  quadratic  fit  to 
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the  first  four  points  is  shown  in  Figure  5.11. 
was  determined  for  each  data  point  by 


Weighting 


w 


i 


t 


where  cn  was  taken  to  be  the  estimated  uncertainty.  Figure 
5.12  shows  a  least-squares  logarithmic  fit,  and  Figure  5.13 
shows  a  combination  of  two  linear  least-squares  fits. 
Results  are  summarized  in  Table  IV. 


TABLE  IV 

SUMMARY  OF  FITTING  FUNCTIONS 

FUNCTION 

V 

THIRD  DEGREE  POLYNOMIAL 

10.894 

QUADRATIC 

8.392 

LOGARITHMIC 

87.632 

DOUBLE  LINEAR  FIRST 

LEG 

3.951 

SECOND 

LEG 

0.749 

O  50  100  150  200  250 


F10.7  (x  1  0“22  Wm“2  Hz'1) 

Figure  5.10.  Weighted  Third  Degree  Polynomial  Fit. 


O  50  100  150  200  250 

F10.7  (xio-22  Wm”a  Hz"1) 

Figure  5.11.  Weighted  Quadratic  Fit. 
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4.  Statement  of  Empirical  Relationship 

Based  on  the  results  of  the  analysis  presented  in 
the  preceding  section,  the  empirical  formula  for  the  scaling 
factor  FRACXR  which  relates  X-ray  flux  and  EUV  flux  in  the 
E-region  of  the  ionosphere  is 

FRACXR  -  14.930  +  0 . 048x( F10 . 7 ) 
for  ri0.7  >  130  and 

FRACXR  -  -24.513  +  0 . 351x( FlO . 7 ) 

for  FlO. 7  <  130,  where  FlO. 7  has  units  10~22  Wm~2  Hz~*. 
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VI .  RESULTS 


In  Chapter  V,  an  empirical  expression  for  FRACXR  was 
proposed.  In  this  chapter,  the  expression  is  incorporated 
into  the  one  dimensional  diffusive  photochemical  model  to 
compute  nitric  oxide  for  two  60  day  time  periods. 

A.  COMPARISON  OF  MODEL  AND  SME  DATA 

Two  60  day  periods  were  chosen  to  further  test  the 

model.  Each  period  was  chosen  to  correspond  to  a  time  where 

the  SME  data  showed  distinct  periodic  variations  in  nitric 

oxide  density.  The  first  60  day  period  was  from  82213  to 

-22  -2 

82273,  and  was  chosen  because  the  F10.7  value  (xlO  W  m 

Hz-1)  ranged  from  a  minimum  of  120  to  a  maximum  of  214  and 

thus  would  correspond  to  FRACXR  values  taken  from  the  upper 

curve  of  Figure  5.13.  The  result  of  this  60  day  run  as  well 

as  the  daily  value  for  F10.7  is  shown  in  Figure  6.1.  Model 

calculations  without  using  X-ray  scaling  are  included  for 

comparison.  The  second  60  day  period  was  from  84210  to 

-22  -2 

84180  ,  and  was  chosen  because  the  F10.7  value  (xlO  W  m 
H3-1)  ranged  from  a  minimum  of  90  to  a  maximum  of  183  and 
thus  would  correspond  to  FRACXR  values  taken  from  the  lower 
curve  of  Figure  5.13.  The  result  of  this  60  day  run,  along 
with  the  F10.7  values  and  unsealed  model  calculations,  is 
shown  in  Figure  6.2. 
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1.  Discussion  o£  Results 


The  average  daily  deviation  of  the  model  calculat¬ 
ions  for  nitric  oxide  density  from  the  SME  satellite 
measurements  was  13.92%  for  the  first  60  day  run  and  24.27% 
for  the  second  60  day  run.  This  is  not  surprising,  since 

the  first  run  was  based  on  a  FRACXR  relationship  which  has  a 

2 

38%  probability  of  being  the  correct  function  (x  *  .75), 

while  the  second  run  was  based  primarily  on  a  FRACXR 

relationship  which  only  has  a  4.8%  probability  of  being  the 

2 

correct  function  (X  *  3.95).  Figures  6.1  and  6.2  show  that 
variations  of  X-ray  flux  by  a  factor  of  about  30  over  the 
solar  cycle  can  explain  the  observed  variation  in  nitric 
oxide  density. 

Calculated  nitric  oxide  densities  for  both  60  day 
runs  are  plotted  against  F10.7  to  test  the  correlation 
between  the  two  data  sets  in  Figure  6.3.  An  inspection  of 
the  figure  shows  that  correlation  does  exist  and  follows  the 
same  "double  linear"  type  of  relation  as  exists  between 
FRACXR  and  F10.7.  A  least-squares  fit  to  each  branch  of 
this  scatter  plot  yields  lines  with  regression  coefficients 
of  0.55  and  0.90  for  the  upper  (F10.7  greater  than  130)  and 
lower  (F10.7  less  than  130)  parts,  respectively.  This 
strongly  suggests  that  FRACXR  is  incorrectly  calculated  for 
values  of  F10.7  greater  than  130. 

An  inspection  of  Figures  6.1  and  6.2  reveals  that 
the  derived  relationship  between  FRACXR  and  F10.7  does  not 


NO  Density  (xIO7 


Figure  6.3.  Correlation  of  Calculated  NO 
Densities  and  F10.7. 


the  model  to  correctly  calculate  NO  densities  for  values  of 

F10.7  ( xIO-22  W  m-2  Hz-1)  from  about  140  to  160.  The 

computed  NO  densities  for  these  values  are  too  high  due  to 

2 

the  artificially  high  production  rate  of  N(  D)  from  the 
increased  X-ray  flux.  This  suggests  that  a  smoother 
transition  exists  between  the  lower  and  upper  curves  of 
Figure  5.13  than  proposed  by  the  "double  linear"  fit. 
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VII.  CONCLUSION 


A.  SUMMARY 

A  theoretical  and  observational  study  of  the  production 
of  nitric  oxide  in  the  E-region  of  the  ionosphere  has  been 
presented.  The  observational  basis  was  the  NO  measurements 
made  by  the  Solar  Mesosphere  Explorer  satellite.  An 
empirical  expression  for  the  variation  in  soft  solar  x-ray 
flux  as  a  function  of  the  10.7  cm  solar  flux  was  determined. 
This  variation  was  in  the  form  of  a  multiplicative  scaling 
factor,  FRACXR ,  between  X-ray  and  EUV  flux  given  by 

FRACXR  *=  14.930  +  0 . 048x(  FlO  .  7  ) 
for  FlO. 7  greater  than  130  and 

FRACXR  -  -24.513  +  0 . 351x( FlO . 7 ) 

for  FlO. 7  less  than  or  equal  to  130,  where  FlO. 7  has  units 
-22  -2  -1 

10  Wm  Hz  .  The  10.7  cm  flux  was  chosen  as  the 
independent  parameter  because  it  is  used  in  the  neutral 
atmosphere  model  as  an  indicator  of  solar  EUV  radiation. 

The  expression  which  calculates  FRACXR  was  incorporated 
into  a  one  dimensional  diffusive  photochemical  model  to 
compute  nitric  oxide  densities.  Model  calculations  for  NO 
densities  showed  significant  improvement  in  both  magnitude 
and  variation.  These  model  calculations  were  compared  to 
observational  measurements  from  the  SME  satellite.  Results 
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of  this  comparison  showed  that  model  calculations  using  the 
above  X-ray  scaling  factor  were  in  much  better  agreement 
with  the  observed  SME  data  than  the  calculations  made 
without  X-ray  scaling.  The  conclusion  from  this  invest¬ 
igation  is  that  variations  of  soft  solar  X-ray  flux  by  a 
factor  of  30  over  the  solar  cycle  can  explain  the  observed 
variation  in  nitric  oxide  densities. 

B.  SUGGESTIONS  FOR  FURTHER  RESEARCH 

Figure  7.1  shows  a  scatter  plot  of  observed  nitric  oxide 
densities  versus  10.7  cm  solar  flux.  Barth  et  al.  (1988) 
pointed  out  that  the  observed  nitric  oxide  densities  show  a 
strong  correlation  with  Lyman-a  flux  over  solar  cycle.  This 
is  shown  in  Figure  7.2.  Examination  of  these  two  figures 
suggests  a  stronger  correlation  of  NO  density  to  Lyman-a 
than  to  F10.7.  It  is  suggested  that  further  investigations 
be  made  into  model  calculations  of  NO  density  based  on  an 
X-ray  scaling  factor  which  is  a  function  of  Lyman-a. 
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Figure  7.1.  Correlation  of  NO  and  F10.7. 


2  2.5  3  3.5  4  4.5 

Lyman— a  (xIO11  photons  cm-2  sec-1) 

Figure  7.2.  Correlation  of  NO  and  Lyman-a. 
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APPENDIX 


NITRIC  OXIDE  REACTIONS  IN  THE  LOWER  THERMOSPHERE 


This  appendix  lists  all  of  the  reactions  used  in  the 
Cleary  one  dimensional  photochemical  model  (1985). 
Reactions  between  neutral  species  are  listed  first,  followed 
by  ion-neutral  reactions.  Equation  numbers  shown  here 
reflect  the  scheme  used  in  Cleary  (1986),  except  for  R39 


which  was  added  later. 

A.  NEUTRAL  SPECIES 

1.  Photoionization 
r 

0  +  hv  ■»  0+  +  e  Rl 

,  o2  +  hv  »►  02+  +  e  R2 

N2  +  hv  N2+  +  e  R3 

O  +  hv  *  0(2D)+  +  e  R4 

N(  4S )  +  hv  N+  +  e  R8 

02  +  hv  ■+  0  +  0+  +  e  R9 

N2  +  hv  -♦  N(4S)  +  N+  +  e  RIO 

NO  +  hv  *►  NO+  +  e  Rl  3 

2.  Photodissociation 

N2  +  hv  *  N(4S)  +  N(4S)  Rl 4 

N2  +  hv  «*  N(2D)  +  N(2D)  R15 

r  NO  +  hv  ■»  O  +  N(4S)  R16 
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3.  Radiative  Relaxation 


N(2D)  *  N(4S)  +  hv 

R32 

4. 

Quenching 

N(2D)  +  e  *  N( 4S)  +  e 

R30 

N(2D)  +  0  •»  N(4S)  +  0 

R31 

5. 

Bimolecular  Reaction 

N{4S)  +  02  NO  +  0 

R28 

N(2D)  +  02  ■»  NO  +  0 

R29 

N(2P)  +  02  NO  +  0 

R39 

N(4S)  +  NO  ■*  N2  +  0 

R33 

N(2D)  +  NO  «►  N2  +  0 

R34 

6. 

Energy  Transfer  or  Photoelectron  Impact 

★  + 

0  +  e  -*0  +  2e 

R5 

02  +  e*  •*  02+  +  2e 

R6 

N2  +  e*  N2+  +  2e 

R7 

N2  +  e*  +  N( 4S)  +  N( 4S)  +  e 

Rll 

N2  +  e*  +  N(2D)  +  N{2D)  +  e 

R12 

IONS 

1. 

Charge  Transfer 

N2+  +  02  *  02+  +  N2 

R18 

0+  +  02  02+  +  0 

R21 

02+  +  NO  NO*  +  02 

R22 

0(2D)+  +  N2  «►  N2+  +  0 

R35 

0(2D)+  +  o2  *  02+  +  0 

R36 
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2.  Ion-Atom  Interchange 

N2+  +  0  +  N0+  +  N(2D)  R17 
0+  +  N2  NO+  +  N(4S)  R20 
02+  +  N(4S)  *►  NO+  +  0  R23 
N+  +  02  •*  N0+  +  0  R25 

3.  Dissociative  Recombination 

N2+  +  e  N(2D)  +  N(2D)  R19 
02+  +  e  +  0  +  0  R24 
NO+  +  e  ■»  O  +  N(4S)  R26 
NO+  +  e  *  0  +  N(2D)  R27 
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